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ABSTRACT 

Using theoretical formalism of F. Dalfovo et. al. condensate 
fraction of BEC (N(/N) as a function of (T /T c °) were evaluated 
using non-interacting and interacting models for two different 
traps and compared with scaling parameter r|. Our results indicate 
that two very different configuration can give rise to the same 
thermodynamic behavior of the condensate corresponding to the 
scaling parameter r|. 

Keywords: Bose-Einstein Condensation, Condensate, Scaling 
Parameter, Thermodynamic Behavior, Interacting Models. 



1. INTRODUCTION 

BEC was observed in remarkable 
series of experiments on vapor of rubidium 
sodium 2 and littium. 3 In all these 
experiments atoms were confined in 
magnetic traps and cooled down to 
extremely low temperature, of the order of 
fraction of micro Kelvin. The first evidence 
for condensation emerged from time of 



flight measurements. The atoms were left to 
expand- by switching off the confining trap 
and then imaged with optical method. A 
sharp peak in the velocity distribution was 
then observed below a certain critical 
temperature providing a clear signature of 
BEC. 

One of -the most relevant features of 
these trapped Bose gases is that they are 
inhomogeneous and finite sized systems. 
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The number of atom ranges typically from 
few thousands to several million. In most 
cases, the confing traps are well 
approximated by harmonic potential. The 
trapping frequency 0)h o provides a 
characteristic length scale of the system a ho 



m(O h 



1/2 



which is in the order of few 
microns. There is a major difference 
between BEC of uniform Bose gas (like 
super fluid helium) and BEC of trapped 
gases. In the case of super fluid helium the 
effect of in homogeneity takes place on a 
microscopic scale fixed by the inter-atomic 
distance whereas in the case of trapped Bose 
gas the size of the system is enlarged as an 
effect of repulsive two body forces. 
Therefore the trapped gases can become 
almost macroscopic objects, directly 
measureable with optical methods. The fact 
that these gases are highly inhomogeneous 
has several important consequences. First 
BEC shown up not in the momentum space, 
as happens in super fluid helium,- but also in 
coordinate space. This double possibility of 
investigating the effect of condensation is 
very interesting from both the theoretical 
and experimental point of view. This 
provides novel methods of investigation for 
relevant quantities like the temperature 
dependence of the condensate, energy and 
density distribution, interference 
phenomenon, frequencies of collective 
excitations and so on. 4 

2. MATHEMATICAL FORMULAE 
USED IN THE EVALUATION 

The Condensate of non-interacting bosons 

One starts with taking the confining 
potential in a quadratic form 



V e Jr) = ™(CD 2 x X 2 +CD 2 y y 2 +C0 2 z Z 2 ) 



(1) 



Here one neglects the atom-atom 
interactions. The many-body Hamiltonian is 
the sum of single-particle Hamiltonians 
whose Eigen values have the form 



(2) 



Where {n x ,n y ,n z } are non-negative integers. 
The ground state (]) (ri...r N ) of N non- 
interacting bosons confined by the potential 
(1) is obtained by putting all the particles in 
the lowest single-particle state {n x = n y = n z 
= 0), namely (j)(ri...r N ) = Yli % where cp 0 (r) is 
given by 



exp 



m 
2h 



(Q} x x 2 + Q} y y 2 +0) z z 2 ) 



(3) 



One has introduced the geometric average of 
the oscillator frequencies: 



a> ho =(a> x a> y a> t y 3 



(4) 



The density distribution then becomes 
n( r ) = N | | <t> 0 ( r ) | 2 and its value grows 
with N. The size of the cloud is independent 
of N and is fixed by the harmonic oscillator 
length : 



a ho = 



y/ 2 



(5) 
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Which corresponds to the average width of 
the Gaussian in Eq. (3). 

Trapped bosons at finite temperature 
(thermodynamic limit) 



In order to evaluate this sum explicitly, one 
usually assumes that the level spacing 
becomes smaller and smaller when N— >oo, 
so that the sum can be replaced by an 
integral 



At temperature T the total number of 
particle is given, in the grand Canonical 
ensemble, by the sum 



n = z H4* w -*<)]- if 



(6) 



While the total energy is given by 

E = „,.{ex P [^K ,,,,-//)]- l} 



(7) 



Where u is the chemical potential and p = 
(KpT) '. Below a given temperature the 
population of the lowest state becomes 
macroscopic and this corresponds to the 
onset of Bose-Einstein condensation. 

The chemical potential is given by 



(8) 



u^> ju, =—haj 
2 

Where 

#7 = (co r + CO v +CO )/3 

v y z > is the arithmetic 

average of the trapping frequencies. 

Inserting this value in the rest of the of the 

sum, one can write 



N-N Q = X 



1 



n x n y n z ±o CXpi 



CO n +CQ n + CO 



•OH 

(9) 



N 



~ N o=f 

Jo 



dndn.dn 



exp: 



+ C0 y n y +C0 



-.)]-> 

(10) 



This assumption corresponds to a semi- 
classical description of the excited states. Its 
validity implies that the relevant excitation 
energies, contributing to the sum (9) are 
much larger than the level spacing fixed by 
the oscillator frequencies. The accuracy of 
the semi -classical approximation (10) is 
expected to be good if the number of trapped 

atoms is large and KpT» ^ 0)h o 

The integral (10) can be easily 
calculated by changing variables 

(/^^=^-). one finds 



K B T 



V 



tico 



N-N 0 ={{3) 

! 7 // t 

(11) 

Where ^(n) is the Riemann ^ function and 
(Oho is the geometric average (4). From this 
result one can also obtain the transition 
temperature for Bose-Einstein condensation 
by imposing that N 0 — > 0 at the transition, 
one gets 



k B T = hO) ho 



( N V " 



<r(3) 



= 0.94fi co ho N 1/3 



(12) 



For temperatures higher than T c ° into the 
chemical potential is less than |i c and 
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becomes N dependent, while the population 
of the lowest state is the order of 1 instead of 
N. The proper thermodynamic limit for these 
systems is obtained by letting N— >oo and co ho 
— ><x>, while keeping the product Noo 3 ho 
constant. With this definition the transition 
temperature (12) is well defined in the 
thermodynamic limit. Inserting the above 
expression for T c ° into Eq. (11) and (12) one 
gets the T dependence of the condensate 
fraction for T< T° 



Ho 

N 



= 1- 



( j V 



\*c J 



(13) 



The same result can be also obtained 
by rewriting Eq. (10) as an integral over the 
energy, in the form 



N 



Jo 



p(e)d£ 
exp((5e) - 1 



(14) 



Putting density of state 
1 

2 



p(e) = j : (ha) ho r 3 £ 2 



From eq n (14), we have 
1 

-ri 



N-N n =f 

Jo 



(hco ho ) 3 e 2 de 



e*-l 



1 1 r £ 2 d£ 



2 (hco 



r t at 

ho) 0^ 1 



1 



1 



We have 



y 2 dy 



N-N n 



=r(3)^(3) = 2^(3) 



hco 



2<f(3) 



ho J 



r K„T^ 



\h®hoJ 



which is the same as that of equation (11) 

On plotting Condensate function (No/N) as a 
function of (T/ T c °) it was found that 
agreement with experimental result 5 is very 
poor. Then two corrections were made in 
order to match the result with the 
experiment. One is finite size effect and 
other is role of dimensionality. 

The finite size correction to the 
equation (13) for the condensate fraction can 
be evaluated analytically by studying the 
large N-limit 6 " 13 . The result for (N 0 (T)/N) is 
given by 



N 



= 1- 



2/3 



AT-V3 
(15) 



The result is shown in table Tj. To 
the lowest order, finite Size effect decrease 
as N 1 13 and depend on the ratio of the 
arithmetic mean 03 and geometric (0) ho ) 
average of the oscillator frequency. For 
axially symmetric traps the ratio depends on 
the deformation parameter 

A = eo z /cD t asdJ/cD ho = (A + 2)/(3A y3 ) 
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Interacting model 

As we know that the effects of two 
body interactions in dilute Bose gas is 
expected to be significant in the presence of 
condensate. In this presence of attraction 
forces there is narrowing of the peak in the 
center of the trap and also increase of peak 
density, one namely discuss the system 
composed of large number of particles 
interacting "with repulsive forces. 

In the presence of repulsive force 
the chemical potential provides an important 
scale of energy lying between oscillator 
energy and the critical temperature 

%0) ho <jU< KgT c ° 

p . One defines a very 

useful parameter r\ which is the ratio 



drdp 



exp 



K p T 



(17) 



Using the Thomas-Fermi approximation for 
the effective mean field potential, one can 
obtain the result 



14 



^ = 1- 
N 



<f(3) 



(18) 



This equation is valid to lowest order in the 
interaction parameter t|. t is the reduced 
temperature t = (T / T c °). 



a 



P ± c 



N' 



V 



a 



K J 



(16) 



This is ratio between the chemical potential 
calculated at T = 0 in TF approximation and 
the critical temperature for non-interacting 
particle in the same trap. 

2 

is a numerical coefficient (N a/a ho ) accounts 
for the effects of two body interactions. If 
one uses the typical values of the parameters 
of the current experiments, One finds that r\ 
ranges from 0.35 to 0.40. Thus, one expects 

that there is an interaction effects of e . 
Using the semi-classical results 14 , one can 
calculate the number of atoms out of the 
condensate. 



3. DISCUSSION OF RESULTS 

In this paper, we have evaluated the 
condensate fraction (No/ N) as a function of 
(T / T c °) using the theoretical formalism 
developed by F Dalfovo etal 14 . In table T : 
and T 2 we have shown the condensate 
fraction (No/N) for interacting and non- 
interacting particles in the different traps. 
For interacting model, we have taken the 
self-consistent mean field calculation within 
Popov approximation for N = 5xl0 4 
rubidium atoms in a trap with a/a ho =5xl0 3 
and X =V8. For N=5xl0 7 Sodium atoms in a 
trap with the value of a/a ho = 5x1 0~ 3 and X 
=0.05. The numerical results were com 
pared with the prediction of scaling theory 
for r| = 0.4. The above calculation indicates 
(in an explicit way) how very different 
configuration can given rise to the same 
thermodynamic behavior if the 
corresponding scaling parameter n is the 
same. It is also interesting to notice that the 



Journal of Pure Applied and Industrial Physics Vol.2, Issue 3, 1 July, 2012, Pages (142-285) 



217 Tejendra Kumar, et al, J. Pure Appl. & Ind. Phys. Vol.2 (3), 212-220 (2012) 



scaling behavior is reached faster in the 
presence of two body interactions than for 
non-interacting particles. 

In table T 3 , we have given the 
evaluated result of the condensate fraction as 
a function of (T / T c °) for axially symmetric 
trap using equation (15) and then the 
obtained result is compared with eq n (13). 
The condensate fraction obtained in this 
way, turns out to be smaller than the 
thermodynamic limit prediction of eq n (13). 
It appears that finite N makes the system 
potentially richer, because new interesting 
regimes can be explored even in case where 
there is no real phase transition in the 
thermodynamic regime. With the numerical 
calculation Ketterle and Van Druten 11 found 
that finite size effects are significant only for 
small values of N less than about 10 4 . 

In table T 4 , we have shown the 
evaluated result of condensate fraction 



(No/N) as a function of t (= T / T c °) for 
various values of 11 (a) n = 0.4, (b) r| = 0.6 
and (c) n = 0 (non-interacting). Result were 
compared with (d) Monte Carlo Calculation 
by Kranth 16 

and experimental data by Ensher et.al 5 . The 
Monte Carlo simulation has been done with 
10,000 atoms interacting through an hard 
core potential. Monte Carlo result 
corresponds to n = 0.35. The experimental 
data corresponds to the n ranges from 0.39 - 
0.45. Monte Carlo result is very much near 
to mean field prediction of n equal to 0.4. A 
detailed comparison between Monte Carlo 
results and mean filed theory have 
performed by Holzmann, Krauth and 
Maraschewski including the analysis of the 
density profiles of the gas at different 
temperatures . Recent calculations also 
confirms the same behavior. 



Table T x 



Evaluated result of condensate fraction (N 0 / N) as a function of (T / T c ) for interacting and 
non-interacting particles for Rb atoms N = 5xl0 4 , a/a h0 = 5.4 xlO" 3 1 = V8, u. = 0.4 



(T/T c °) 


(No/N) 


Interacting 


Non-Interacting 


0.1 


0.787 






0.968 


0.2 


0.924 






0.942 


0.25 


0.896 






0.905 


0.30 


0.824 






0.832 


0.35 


0.745 






0.738 


0.40 


0.692 






0.674 


0.45 


0.612 






0.625 


0.50 


0.554 






0.567 


0.55 


0.526 






0.567 


0.60 


0.497 






0.506 


0.65 


0.432 






0.414 


0.70 


0.386 






0.408 


0.75 


0.339 






0.367 


0.80 


0.268 






0.318 
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Table T 2 



Evaluated result of condensate fraction (N 0 / N) as a function of (T / T c °) for interacting and 
non-interacting particles for Sodium Atom N = 5xl0 7 a/a h0 = 1-7 xlO" 3 , k = 0.05, \\ = 0.4 



(T/T c °) 


(No/N) 


Interacting 


Non-Interacting 


0.1 


0.992 






0.982 


0.2 


0.956 






0.946 


0.3 


0.897 






0.875 


0.4 


0.824 






0.806 


0.5 


0.786 






0.765 


0.6 


0.697 






0.669 


0.7 


0.547 






0.525 


0.8 


0.512 






0.452 


0.9 


0.467 






0.408 


0.95 


0.418 






0.382 


1.0 


0.359 






0.339 



Table T 3 



Evaluated result of condensate fraction (N 0 / N) as a function of (T / T c ) for axially 
symmetric traps using equation (15). The results were compared with equation (13) 





(No/N) 


(T/T c °) 


Using eq° (15) taking finite size 


Using eq n (13) without finite 




effect 


size effect 


0.1 


0.967 


0.982 


0.2 


0.942 


0.956 


0.3 


0.918 


0.942 


0.4 


0.894 


0.908 


0.45 


0.813 


0.832 


0.50 


0.753 


0.765 


0.55 


0.719 


0.735 


0.60 


0.627 


0.644 


0.65 


0.584 


0.622 


0.70 


0.548 


0.608 


0.75 


0.502 


0.554 


0.80 


0.343 


0.468 


0.85 


0.309 


0.365 


0.90 


0.228 


0.242 


0.95 


0.182 


0.204 
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Table T 4 



Evaluated result of condensate fraction (No/N) as a function of (T / T c ) for various values of 
x\ (a) Ti = 0, (b) i\ = 0.4, (c) i] = 0.6 and (d) Monte Carlo result (e) Experimental result. 





(No/N) 


(T/T c °) 


n = 0 


r\ = 0.4 


n = 0.6 


Monte Carlo 
results 


Expt. 
Results 


0.1 


0.995 


0.942 


0.922 


0.953 


0.939 


0.2 


0.976 


0.910 


0.894 


0.922 


0.916 


0.3 


0.942 


0.886 


0.856 


0.896 


0.890 


0.4 


0.896 


0.823 


0.796 


0.834 


0.842 


0.45 


0.853 


0.785 


0.763 


0.792 


0.788 


0.50 


0.832 


0.712 


0.709 


0.724 


0.709 


0.55 


0.816 


0.695 


0.675 


0.708 


0.686 


0.60 


0.802 


0.638 


0.612 


0.644 


0.642 


0.65 


0.786 


0.602 


0.588 


0.616 


0.619 


0.70 


0.769 


0.556 


0.543 


0.569 


0.573 


0.75 


0.705 


0.514 


0.496 


0.526 


0.532 


0.80 


0.682 


0.448 


0.428 


0.455 


0.431 


0.85 


0.605 


0.366 


0.343 


0.375 


0.374 


0.90 


0.538 


0.298 


0.277 


0.302 


0.308 


0.95 


0.233 


0.124 


0.116 


0.133 


0.122 
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